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Abstract 

Focused ion beam (FIB) tomography provides high resolution volumet¬ 
ric images on a micro scale. However, due to the physical acquisition 
process the resulting images are often corrupted by a so-called curtain¬ 
ing or waterfall effect. In this paper, a new convex variational model 
for removing such effects is proposed. More precisely, an infimal convo¬ 
lution model is applied to split the corrupted 3D image into the clean 
image and two types of corruptions, namely a striped part and a lam¬ 
inar one. As regularizing terms different direction dependent first and 
second order differences are used to cope with the specific structure 
of the corruptions. This generalizes discrete unidirectional total vari¬ 
ational (TV) approaches. A minimizer of the model is computed by 
well-known primal dual techniques. Numerical examples show the very 
good performance of our new method for artificial and real-world data. 
Besides FIB tomography, we have also successfully applied our tech¬ 
nique for the removal of pure stripes in Moderate Resolution Imaging 
Spectroradiometer (MODIS) data. 


1 Introduction 

The motivation for the present work was the analysis of aluminum matrix 
composites reinforced with silicon carbide particles by given high resolution 
3D FIB tomography images. 

FIB tomography, also known as serial sectioning, is an imaging technique, 
that has been used commercially for about 20 years. It is applied for prepa¬ 
ration and direct observation of structural cross-sections and the generation 
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Figure 1: Illustration of a FIB system. A focused ion beam mills slice by 
slice of the sample in z-direction and the exposed surface is imaged by an 
electron beam. Image credit: [H] 


of microstructural data in three dimensions. While classical A-ray tomogra¬ 
phy does often not reach the required material resolution, FIB tomography 
enables to investigate structures on a scale down to several nanometers. 
The FIB system combines a classical scanning electron microscope with 
an ion beam. An illustration is shown in Fig. The ion beam mills and 
polishes the material sectionally, while the electron beam is used for imaging 
the surface after each section. Several hundred of these serial slices finally 
form a 3D image. For a more detailed introduction to FIB tomography we 
refer to [a uni [29]. 

Unfortunately, such images often suffer from the so-called curtaining effect 
HHnsniisi, see Fig. This effect arises because the sputtering rates of the 
ion beam are sensitive to local changes in the surface structure. In particular 
below pores or cracks the rates vary and lead to stripe-like artifacts in the 
images. For the aluminum matrix composite in Fig. e.g., at the phase 
boundary between aluminum and the reinforced silicon carbide particles, due 
to their locally varying material characteristics. Further, it may happen that 
the material is milled incompletely. This leads to bright laminar artifacts. 
To reduce the curtaining effects, a lower beam current could be used but this 
leads to a considerably higher milling time. Note that the whole imaging 
and milling process already takes about 25 hours for each data set considered 
here. So further reducing the beam current leads to impractical processing 
times. 

Therefore, we propose a computational model to reduce curtaining effects. 
The corruptions in our data consist both of stripe-like structures due to 
different sputtering rates and laminar structures due to incomplete milling 
of the material. The corrupted parts cannot be used for the further analysis 
of the material unless the curtaining effects are removed. For our specimen 
in Fig.j^this was the case for more than half of the data. In FIB tomographic 
images the laminar corruptions are often cut off as they only occur in the 
lower part of the image. However, we will keep them since it is not desirable 
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Figure 2: Volume image of an aluminum matrix composite obtained by FIB 
tomography (255 x 255 x 100 pixels, i.e., 100 slices of size ll.G/rmx 11.6//m) 
and typical slices in the x-y, x-z and y-z plane. 

to throw away half of the data. 

The problem of removing stripes in images was tackled by various ap¬ 
proaches: Fourier based filtering was suggested in [ninij, moment matching 
in [18] and histogram based methods in [241135] . Recently, variational mod¬ 
els were successfully applied for the destriping of 2D images. The central 
idea in |5| is the application of unidirectional discrete TV terms to extract 
the stripes by minimizing the functional 

argmin ||Vy(/ - u)||i -k A||V 2 ;u||i, A > 0, 

U 

where / is the original image and u the destriped one. This method was 
improved in [T2| by adding a least squares data term and a framelet regu¬ 
larization term: 

argmin ^11/ - u\\l + i^i\\Vy{f - u)||i -k i' 2 \\^xu\\i + i/ 3 ||VFn||i, (1) 

where W denotes a Parseval framelet transform. It has been shown that this 
model gives much better results than the previous one also without framelet 
regularization, i.e. t '3 = 0. The framelet regularization led to another slight 
improvement. Finally, a more general method for the variational denoising 
of images with structured noise was developed in which has also been 
applied for destriping. 

Applications of the aforementioned destriping methods are, e.g., the restora¬ 
tion of MODIS data jS] IT^ [35] and FIB tomography images [IlKIZj. In this 
paper, we focus mainly on the latter, where our FIB images are not only 
corrupted by stripe structured noise, but also by laminar effects. To the 
best of our knowledge the removal of curtaining effects involving additional 
large laminar parts was not considered before. 

Compared to the above destriping methods our new model takes the follow¬ 
ing aspects into account: 
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Al) The existing destriping methods remove stripes only. Since our data 
is not only corrupted by stripes, we split the corrupted image not only 
into a clean part and stripes, but also into a third, laminar part. 

A2) To avoid smoothing and to keep also the small image details, we use 
a hard constraint for the decomposition of the corrupted image into 
the clean, the striped and the laminar part by applying an infimal 
convolution model. 

A3) We consider 3D data whereas most destriping methods, except [T7] . 
focus on 2D images. 

To deal with A2) we propose an infimal convolution model for the splitting. 
The infimal convolution was applied for the first time in image processing 
by Chambolle and Lions followed by many other papers, see, e.g., the 
general discrete approach in 133 or the continuous function space approach 
in [53] • We also refer to a nice PhD thesis on infimal convolutions m- 
Infimal convolution models were in particular used for the decomposition of 
images. 

Various variational decomposition models were proposed in the literature. 
The additive splitting into geometric and oscillating parts as texture and 
white Gaussian noise typically exploits TV- or Besov-seminorms for the first 
one and Meyers’s G-norm jT] |28l |39|, the norm of the dual Sobolev space 
|3QII32]) the squared L 2 -norm of the DCT-transformed texture [38] 
for the second one. For specifying the noise component also the norm of the 
dual Besov space ^ was successfully used in [2] . Decomposition models 
were also applied for simultaneous structure-texture inpainting, e.g., in [3] 
as an extension of the so-called morphological component analysis model 
[38] . For the separation of point and curved structures based on wavelets 
(Besov norms) and shearlets or curvelets we refer to [261138j . The clou of 
all these methods is the adaptation of the additive components to the task 
at hand. For the curtaining problem we will show that our directional TV 
model is very well suited. 

Organization of the Paper. In Section [^ we introduce our variational 
model for the removal of curtaining effects and motivate its choice by the 
special structure of the corruptions. Then, in Section [^ a primal dual 
algorithm is presented for finding a minimizer that consists of the clean 
image and two types of corruptions. Section|^demonstrates the performance 
of our algorithm for volume images obtained by FIB tomography as well as 
artificial examples. Further, a comparison to existing methods for destriping 
and denoising of MODIS data is given. Finally, conclusions are drawn in 
Section [5] 
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2 New Model 


We are given a 3D image F £ [0, corrupted by curtaining effects. 

The total corruption can be seen as an overlay of two parts: 

i) Laminar part consisting of very bright areas in the x-y-plane. Along 
the z-axis these regions occur very fluctuant and it is very unlikely 
that the effect appears in the same region in two subsequent frames. 

ii) Stripes along the y-axis which are few pixels wide in x and z direction. 

Based on the previous observations the corrupted image F is assumed to be 
an additive composition of the clean image U, the stripes S and the laminar 
part of the corruptions L, i.e., F = U + S + L. 

To give a sound discrete formulation of our model which takes i) and ii) 
into account, all images are considered as column- and slice-wise reshaped 
vectors in with N := UxUyUz- The reshaped images are denoted by small 
letters so that in particular / = u + s + 1. Then all linear operators can 
be written as matrices and their adjoints are just the transposed matrices. 
By (8) we denote the Kronecker product of matrices. We use the following 
directional forward difference operators: 

Vy:=Dy, V 

_ 

^X,y,Z ■ - I Fy I , .- Dzzj 

\DzJ 

where 

Fx ■ — Iriz ® ^riy ® Friz ;) Fy ■ — In^ ® Friy ® Inx i 

Fz •= Friz ® ^n-y ® Irix i Fzz ■— F^^ ® luy ® Irir 

and 

\ 

e 

1 

0 / 


-2 1 
0 0 

Note that we assume constant extension of the image beyond the bound¬ 
ary for the first order differences and linear extension for the second order 
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differences which is exactly described by the matrix form of the operators. 
Let us further emphasize that we use the Kronecker product notation only 
for a correct description of our model in a convenient matrix-vector form. 
In our numerical computations we work with arrays based on the relation 

\ec{AXB'^) = {B ® A) vec(X), 

for A G B G F G 

For a vector w G appearing, e.g., with d = 2 after application of Vx,z 
or Vx,y and d = 3 after applying Vx,y,zi the grouped £ 2 ( 1 ^'^) — norm 

is given by 


N 

:= ^ 




d-l 

E 

j=0 


w. 


i+jN' 


Finally, let lq denote the indicator function of a set C defined by 

^ ^ _ fo if tt G C, 

1 -boo otherwise. 

We propose the following infimal convolution model 

inf ipi{u)Fip 2 {s) + ip^{l), 

u-\-s-\-l=j 


( 2 ) 


with 


ipi{u) := /Ui||Va;, 2 u|| 2 ,l + lX 2 \\^zU\\l + i[o,l]Jv(u), 

^ 2 {s) := ||Vj^s||i, 

and regularization parameters //i,^ 2 ,I ^3 > 0. For convenience, we also add 
the definition of the above norms for the corresponding array formulation; 

>^2 1 

\\'^x,zU\\2,i '■= ^ {Ui+i^j^k — Uij^k)"^ + {Uij^k+1 — , 

rix ,TLy ,nz 

IIA^t/lli := ^ \Uij^k+i — 2Uij^k + Uij^k-i\, 

Tlx ,TLy ,nz 

l!VyS'||i:= \Uij+i^k — Uij^kl 

ij,k—l 

with the appropriate mirrored or linear extensions at the boundary. The 
last summand in (pi ensures that the clean image u has its range in [0,1]. 
The choice of the other terms is motivated by i) and ii) as follows: 
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i) II Va;,y^|| 2 ,i and HA^nUi: laminar distortions have a small width in z di¬ 
rection, thus a high response in the second order derivative is expected 
along this direction. Therefore, we penalize second order differences 
in z direction of u. Additionally, these corruptions are spacious in x 
and y direction such that a bidirectional TV regularizer is useful for 1. 

ii) ||Vj;s||i and ||Va;, 2 u|| 2 ,i: stripes occur only in y direction so that we 
penalize first order differences in y direction of s. This was also done 
in mm setting s = f — u. Further, we penalize coupled first order 
differences in x and z direction of u since we assume that the main 
differences along these directions occur in s while u has only few edges. 

Remark 1. Alternatively to //i||Va;^ 2 i^|| 2 ,i + /^ 2 ||A 2 u||i we can also use the 
3D TV term \\S/x,y,z'U‘\\ 2 ,i, he., instead o/ Q we are asking for 

inf ifi{u) + (p 2 {s) + if 3 {l), (3) 

R+S+/ = / 

where ipi, i = 2,3 are defined as in (§ and 

;= lJLl\\Vx,y,zU\\2,l + t[ 0 ,l]jv(u). 

The term for u corresponds to the assumption that u is a piecewise smooth 
image. Our model Q focuses more on the corruptions. We argue in i) that, 
in contrast to the laminar corruptions, u has a small second order derivative 
in z direction. This viewpoint is also taken in 13 Cl/. We have implemented 
both methods. The difference between the models is best visible in Fig. in 
the numerical section. 

Proposition!^ shows that our models Q and Q have a solution. Its unique¬ 
ness cannot be expected since arbitrary constants can be added to s and 
subtracted from I without changing the objective function. 

Proposition 1. Let f G [0,1]^. Then respectively Q, has a minimizer. 

The proof is given in the appendix. 

3 Algorithm 

In this section, we propose to use a primal dual algorithm to solve the convex, 
but non-smooth problem (©• Problem @ can be handled in a similar 
way. Primal dual algorithms have recently found wide applications in image 
processing and machine learning. The reason is that these algorithms split 
the original task into a sequence of simple computable, proximal mappings. 
This paper cannot address the huge amount of papers on various primal dual 
strategies, which were successfully applied in imaging. For a good overview 
on primal dual methods we refer to da [25]. 
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Recall that the proximal mapping of a proper, convex, lower semi continuous 
function (p-. —)■ (— 00 , + 00 ] is defined by 

prox;, (x) := argmin | (f{y) + - y\\l 

j/gRiV L 

Indeed the proximal mapping at x G exists and is unique. To make 
primal dual algorithms efficient, the proximal mapping should in addition 
be simply computable. For applications of proximal mappings see, e.g., m- 
We want to apply the primal dual algorithm with an extrapolation of the 
dual variable suggested by Chambolle, Pock and Cremers in [ini[33], see 
also |lll 121113^ for extensions. We call this algorithm PDHG (Primal Dual 
Hybrid Gradient) algorithm in accordance with the name of a similar al¬ 
gorithm which was presented without the extrapolation idea and without 
convergence proof in |44) . Later its convergence was examined in [3]. 

To this end, we have to rewrite our model Q. Using the vectors x := 
(u, s, l)'^ and y := and the notation 

C := {(tt, s,l) ■. u + s + l = f, n G [0,1]}, (4) 

Ky) ■= y^ihih,! + M2II2/2II1 + llyalbp + M3II2/4II1 

we ask for 

arg min {tc(x)-|-/i(y)} subject to Kx = y, 

^,y 

where 


(Vx,Z 

0 

0 

\ 

A, 

0 

0 


0 

Vy 

0 


V 0 

0 

^x,y/ 



Then the PDHG algorithm from [101133] with an extrapolation of the dual 
variable p reads 

= prox.,,^(x^ - tK^P^) = Vcix’^ - 
/+i = prox,,*(/ + uiLx^+i), 

/+i=/+i+0(/+i-/), 0G(O,1], 

where h* denotes the conjugate function of h and Vc is the orthogonal 
projection onto C. The algorithm converges if the parameters fulfill ra < 
. For our structured matrix K it is not hard to check by diagonalizing 

it with trigonometric transforms that ||iL ||2 < 24, so that xa < ^ must be 
fulfilled. 
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Applying Moreau’s identity, one can directly work with the function h in¬ 
stead of its conjugate. Using the rescaling := the algorithm becomes 

Vc{x^ - 

proxi^(6*’ + A:x^+^), 

, 

0E(O,1]. 

For details, see, e.g., [8]. Fortunately, the proximal mapping of h can be 
separated into the proximal mappings of its summands with respect to yi, 
i = 1,..., 4. The whole algorithm with our parameter choice is detailed in 
Algorithm 


yk+l ^ 
6^+1 = 
1^+^ = 


Algorithm 1: PDHG for Q. 

Initialization: = f, = 0, = 0, 6® = 0, = 0, 

Z = 1,2,3,4, 0 = 1, r = |, cr = g. 

Iteration: For k = 0,1,... iterate 

{u,s,lY'^+^^ = - rcjAj6f\ 

= prox^||.||^^^(6f^-b 

^(fc+i) ^ proxA^||.||^(6^^^ -b 
= proxi||.||^(6f^ -b 
= prox^||.||^_^(6f^-b 

^(fc+1) ^ ^{k) ^ ^^^^^{k+1) _ y{k+l) 

6^^+^) = 

6^+1) = 6^ + - yf 

^(k+i) ^ ^(fc) ^ y^^(fc+l) _ y{k+i) 

^(fc+i) ^ ^(fc+i) ^ Q^iy{k+i) _ ^(fc)^ i = 1,..., 4 

Output: Clean image u, corruptions s, I 


The first step of the algorithm is a voxelwise projection Vc onto the set 
C. Since C is just a plane with an additional range constraint for u, this 
projection can be computed in a straightforward way. The update steps for 
y 2 and ya require a componentwise soft shrinkage, which is defined for t G M 
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by 


Sx{t) := 



if |i| < A, 

1^) if 1^1 > A. 


For yi and a componentwise coupled shrinkage procedure is required. For 
vectors t = (ti,t 2 )^ G 1^^ it is given by 


SA(t) 


0 if ||t||2 < A, 

if||f||2>A. 


Due to the separation of the computation of the i = 1,..., 4, the algo¬ 
rithm can be implemented in parallel on a multi-core architecture. Recently 
such parallel strategies were also examined by Pesquet and co-authors [ZllMi. 
Further we want to mention the interesting stochastic block coordinate al¬ 
gorithms in [161ES] ■ 


4 Numerical Examples 

In this section, we present the results of our algorithm, which was im¬ 
plemented in MATLAB2015b, for 3D artificial data, 3D FIB tomography 
data and 2D MODIS data. For the 3D examples, videos of the whole vol¬ 
ume images are provided at our website http;//www.mathematik.uni-kl.de/ 
imagepro/members/fitschen/curtainingO/. 

We compare the results of our model ( [2| w ith three other methods Ml, M2 
and M3 and with the modified model Q in Remark Since there seems 
to be no work especially on curtaining effects, the proposed methods are 
compared to the following 3D destriping and denoising algorithms: 

Ml) Firstly, we generalize the 2D destriping approach from [12] with regu¬ 
larization terms as in 0 without the framelet part to 3D data, which 
results in 

argmin^ll/ - u\\l + I'lWVyif - u)||i -b i^ 2 \\'^x,zu\\ 2 ,i- 

To remove the corruptions that are not stripes, a median filter in z 
direction is applied beforehand since the simple destriping algorithm 
is not able to remove these artifacts. 

M2) Secondly, we consider the 2D denoising method proposed in [T7| for 
removing structured noise. We have used the implementation available 
at http://www.math.univ-toulouse.fr/~weiss/PageCodes.html 
with a Gabor filter and parameters {p,a). Notice that this method 
is applicable for a wide range of structured noise. Thus one can not 
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expect that the results are as good those of our model which is adjusted 
to this specific task. Again, we apply a median filter in 2 direction first 
and then apply the algorithm to the 2D slices separately. 

M3) Thirdly, we consider the 3D shearlet hard thresholding method in [27] 
for denoising of videos, see also [22|. We have used the software pack¬ 
age ’’ShearLab 3D” available at http: //www. shearlab. org/software. 
For this method we have used the shearlet system SL3Di, i.e., 3 scales 
and a redundancy of 76. In total we have used 5 different thresholding 
factors, one factor ui for the direction orthogonal to the laminar part, 
another one V 2 for the directions nearest to that direction, in the same 
way another two for the stripes (^ 3 ,^ 4 ) and finally a hfth factor U 5 for 
all other directions. 

Compared to the methods Ml and M2 our method has the advantage of 
computing everything in one step without any preprocessing. 

Artificial Data 

The first two experiments deal with artificial 3D images disturbed by stripes 
and laminar structures, that are fluctuant in z direction. The parameters 
of all models were chosen by an extensive grid search such that the PSNRs 
were optimized. 

In Fig. a comparison of the proposed models © and i) and the methods 
Ml, M2 and M3 is shown. Specific slices of the volume image are presented. 
The method Q performs clearly better than the other ones. This is also 
confirmed by the corresponding PSNR values for the whole 3D data. If we 
apply Ml and M2 without a median filter in the preprocessing step, the 
PSNR is even worse, namely PSNR = 12.53 for the destriping method and 
PSNR = 10.05 for the denoising method. In contrast to these methods 
with M3 it is possible to remove the laminar part without an additional 
median filter. But nevertheless the results are clearly worse compared to 
the proposed model. 

Comparing the models Q and Q, we observe that the main differences 
occur along the laminar corruptions. To exclude that the reason for this 
behavior are the sharp edges of the laminar part, we show another exam¬ 
ple with a smoothed laminar part in Fig. [^ Model gives an almost 
perfect reconstruction, whereas the method ® is still not able to remove 
the laminar corruptions correctly. The results are summed up in Table 
Besides the PSNR, which we have discussed, the mean square error (MSE) 
and structural similarity index (SSIM) are stated in the table. In all cases 
our method Q performs best and the methods Ml, M2 and M3 are clearly 
worse. 
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corrupted proposed Q proposed ([^ 



Figure 3: Results of various methods for the artificial image. Top left to 
bottom right: Exemplary x-y slice of the corrupted image, reconstructed im¬ 
ages with our method Q ((/ii, |U 2 ,/Us) = (i^, 3 ^, 3 ^), PSNR 34.62), our 
method ^ {{^ 1 ,^. 3 ) = ( 3 ^, 3 §o)> PSNR 27.74), the generalized destriping 
method Ml ((z^i,zz 2 ) = (0.9,0.5), PSNR 14.31), the generalized denoising 
algorithm M2 (p = l,a = 0.13, PSNR 14.05) and the shearlet hard thresh¬ 
olding method (m = 5 ,U 2 = 3 ,U 3 = 1 . 1 , U 4 = 1 ,U 5 = |). 


Image Error measure 




PSNR 

MSE 

SSIM 

PSNR 

MSE 

SSIM 


Model (p|) Model ([3|) 

34.62 27.74 

0 . 35 - 10-3 1 . 68 - 10-3 
0.99 0.959 

32.76 27.39 

0.53-10-3 1.82-10-3 
0.99 0.959 


Ml M2 M3 
14.31 14.05 12.17 

0.0371 0.0394 0.0607 
0.805 0.784 0.741 
14.14 13.92 12.22 

0.0386 0.0406 0.06 
0.808 0.782 0.7512 


Table 1: Error measures for the results of various methods. 
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Figure 4; Results of the proposed methods for an artificial image. Top 
to bottom: exemplary slices in x-y, x-z and y-z directions of the volume 
image. Left to right: original image, corrupted image, reconstructed image 
u with our method Q ((^i,// 2 ,/^s) = 3 ^) 3 ^)) PSNR 32.76) and ^ 

((m,f^3) = (5M>30o)’PSNR 27.39). 
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corrupted reconstructed stripes laminar 



Figure 5: Results of our method Q for real FIB data. Left to right: Exem¬ 
plary x-y slices of the corrupted image F, reconstructed image U, corrup¬ 
tions S and L. The first row shows a slice with prominent corruptions also 
in L and the second row with almost only stripes. 

FIB Data 

Next we apply the methods to different FIB tomography data, namely to 
aluminum matrix composites with different particle sizes and to a bronze 
sample. The images of the aluminum matrix composite were obtained from 
the “Nano Structuring Center” Kaiserslautern and the bronze sample from 
the “Material Engineering Center Saarland” in Saarbriicken. 

For all data sets we applied the proposed model Q with the same parame¬ 
ters, namely 

(/^i) 1^2;= ( 300 ’ 300’ 3Uo)- 

Only for the bronze sample we set ^2 = 0 since there is no laminar part. 
Fig. [^demonstrates that our method Q is able to split the corrupted image 
nicely into a clean image and two corrupted parts. 

In Fig. [^ a comparison of our model Q with Q, Ml and M3 is presented. 
For the slice that is mainly corrupted by stripes, Ml shows acceptable re¬ 
sults, but cannot compete with Q and Q, which perform equally well here. 
If there are not only stripes, one has to apply a median filter beforehand to 
obtain reasonable results with Ml. However, this results in a loss of small 
image structures like the bright aggradations in the upper left part of the 
images. The method M3 leads to slightly smoothed results and still shows 
some remains of both stripes and the laminar part. 

Finally, we show results obtained by our model Q for FIB images of various 
materials. In Fig. exemplary slices are shown. We obtain very good 
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proposed O 


Ml 


M3 



Figure 6 : Results of various methods for real FIB data. Left to right: 
Exemplary x-y slices of the corrupted image, reconstructed image with the 
proposed model Q, the proposed model Q ((/xi,^ 3 ) = ( 3 ^, 3 ^))) the 
destriping method Ml ((z^i,r' 2 ) = (0.4,0.2)) and shearlet thresholding M3 
(t>i = 4,U2 = 3 ,i; 3 = 3,t>4 = 3, Us = ^). 

results although the materials as well as the corruptions look different. For a 
comparison also the results of Ml and M3 are presented. However, these are 
clearly worse. In particular with Ml the laminar structures are still present 
and M3 is not able to remove the stripes. But note that both methods are 
not designed for such corruptions. 

MODIS Data 

with 2D MODIS data (http://ladsweb.nascom. 
also used in [5]) show that our method works also 
The parameters were chosen to give the best visual 
impression. The results of the various methods are presented in Fig. For 
a comparison we have used the full model with the framelet regularization 
Q, proposed by [12], and the method M2. In both examples, all methods 
lead to visually good results. This shows that the proposed method can 
also be applied for destriping, although it is designed for the more general 
setting of stripes in combination with laminar corruptions. 

5 Conclusions 

We proposed a variational model for removing the curtaining effect in FIB 
tomography images. Due to the special structure of the corruptions, we 
chose an infimal convolution model consisting of directional TV terms and 
also second order terms for u. Alternatively, a full 3D TV term for u may be 
applied. The numerical experiments show that the proposed model leads to 


Our final experiments 
nasa.gov/index.html, 
fine for destriping only. 
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Figure 7: Results of our method (§ and Ml {{vi,V 2 ) = (0.4, 0.2)) for differ¬ 
ent materials. From top to bottom: aluminum matrix composite with small 
particles, large particles and a bronze sample. 


very good results. So far, all FIB tomography experiments share the same 
model parameters. This makes the method highly suitable for applications. 
Our restoration method enables us to use the FIB tomography results for 
the further (statistical) analysis of the 3D structure of the material, which 
was up to now only possible for a considerable smaller area of the data. 

For future work we plan to parallelize the algorithm taking also spatial 
domain decompositions into account. The aim is an implementation of our 
approach in a way that it can be run simultaneously to the image acquisition 
process. Here we will also take other algorithms with potential on multicore 
architecture as e.g. [aiMj into account. 

A Appendix 

Proof of Proposition^ We rewrite Q as 

E(u, s, 1) ■.=(pi{u) + ip 2 {s) + (fsil) -h f}{u,S,l). 

A feasible point is obviously given by (u, s, 1) = (/, 0, 0) and E{f, 0, 0) < C 
for some C > 0. 

Assume that the infimum is not attained. Then there exists an inhmal 
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original 


proposed O 



Ml M2 



Figure 8: Results for destriping 2D MODIS data using three different meth¬ 
ods: our proposed method Q with (/ii,/^s) = (0.5,1,4), the destriping 
method Ml [12] in Q with (i^i, 1 / 2 , t's) = (§) |) 1) and the denoising method 
M2 [T7] with p = l,a = 1. 
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sequence {xjjjgis^, where Xi := {ui,Si,li) so that 

E{xi) < E{xi-i) < C, ||a;^||oo —^ oo, E{xi) —)• inf£^(a:) 

X 

as z —>■ oo. We will show that we can find another sequence {xijigN with 
||(5i)^||oo < C, E{xi) = E{xi). 

Since this new sequence is bounded, there exists a converging subsequence. 
By the lower semi-continuity of E this implies that the limit of this subse¬ 
quence is a minimize!’ of E contrary to the assumption. 

It remains to show that there exists Consider an arbitrary element 

X = {u, s, 1) of Let the minimal and maximal value of the j-th x-y- 

plane of I be denoted by mj, respectively Mj. Then we have by E{xi) < C 
and /U 3 > 0 that Mj — rrij < C. Now define 



f mi\ 

/ mi \ 

1 = 1- 

: G) ln,,ny and s = s 



KmnJ 

K'mnJ 


where In^ny = (1, • • • j 1)"^ £ , Obviously, E{u, I, s) = E{u, I, s). Fur¬ 

ther we have ||Z||oo < C and by the constraint on u also ||u||oo < 1- Since 
the given image / is also bounded, we conclude by / = u -|- s -|- / that s is 
bounded, too. This finishes the proof. □ 
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